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Abstract 

Anurans breed in a variety of aquatic habitats with contrasting levels of desicca- 
tion risk, which may result in selection for faster development during larval stages. 
Previous studies suggest that species in ephemeral ponds reduce their develop- 
mental times to minimize desiccation risks, although it is not clear how variation 
in desiccation risk affects developmental strategies in different species. Employing 
a comparative phylogenetic approach including data from published and unpub- 
lished studies encompassing 62 observations across 30 species, we tested if species 
breeding in ephemeral ponds (High risk) develop faster than those from permanent 
ponds (Low risk) and/or show increased developmental plasticity in response to 
drying conditions. Our analyses support shorter developmental times in High risk, 
primarily by decreasing body mass at metamorphosis. Plasticity in developmental 
times was small and did not differ between groups. However, accelerated develop- 
ment in High risk species generally resulted in reduced sizes at metamorphosis, while 
some Low risk species were able compensate this effect by increasing mean growth 
rates. Taken together, our results suggest that plastic responses in species breeding 
in ephemeral ponds are constrained by a general trade-off between development 
and growth rates. 



Introduction 

One of the primary goals of life history theory is to explain the 
variation in age and size of organisms at ontogenetic niche 
transitions (hatching, metamorphosis, maturation) , and how 
potential trade-offs may constrain its plasticity and evolution 
(Wilbur and Collins 1973; Stearns and Koella 1986; Gotthard 
and Nylin 1995). For instance, when the developmental time 
is constrained- for example in amphibian tadpoles developing 
in desiccating ponds (Newman 1992) or insect larvae during 
shortening day lengths (Johansson et al. 2001)-optimality 
models predict that organisms should accelerate their devel- 
opment while increasing growth rates, to minimize juvenile 
mortality and the costs of decreased size at the life-history 
transition (Abrams et al. 1996; Strobbe and Stoks 2004). 
However, a faster development may result in fitness costs 
at other levels, such as decreased survival during adulthood 
(Arendt 1997; Nylin and Gotthard 1998; Stoks et al. 2006), 



hence the adaptive value of this response will depend on the 
balance between its benefits and costs during larval develop- 
ment and the subsequent stages. It remains consequently un- 
clear how species or populations subject to different selection 
intensities for accelerated development eventually respond to 
selection, and how they can minimize the impact of increased 
developmental rates on overall fitness. 

Larval amphibians occupy aquatic environments along a 
wide permanency gradient (Woodward 1983; Richardson 
2001; Babbitt et al. 2003), and they are an ideal system for 
examining the evolution of developmental strategies in re- 
sponse to environmental variability. It is often assumed that 
species breeding in streams and permanent ponds, which are 
subjected to low desiccation risk, have longer developmen- 
tal times than those breeding in temporary and ephemeral 
ponds (Denver 1997; Wells 2007). If tadpoles must develop 
faster to minimize desiccation risk, either as an evolution- 
ary response or by means of developmental plasticity, three 
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different scenarios can emerge (Fig. 1). First, developmental 
period can be modulated by increasing or decreasing body 
mass at metamorphosis without requiring any change in av- 
erage growth rate (this scenario is shown in Fig. la). Second, 
tadpoles can grow faster to compensate for the acceleration 
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Figure 1. Studying the relationship between developmental rates 
(=1/developmental period), mean growth rates, and body mass at 
metamorphosis. Compared to the developmental trajectory under non- 
stressful conditions (i.e., no desiccation) represented by the 1:1 line, 
accelerating development would result in the following scenarios: (a) if 
developmental rates and mean growth rates are uncoupled, for exam- 
ple, growth rates may be constant to differing rates of development and 
thus a decrease in developmental period due to pond drying would be 
accompanied by a similar reduction in size at metamorphosis; (b) if a min- 
imum viable body size is required to enter metamorphosis or for survival 
as an adult, mean growth rates should increase to compensate for the 
accelerated development; and (c) developmental rates would increase 
at the expense of growth if there is a trade-off between these separate 
processes. Note that this simplified scheme depicts a linear growth tra- 
jectory for illustrative purposes because we work with average growth 
rates throughout the study, real growth trajectories in anuran tadpoles 
are by no means linear (see Fig. S1). 



in development, thus maximizing size at metamorphosis. 
Such a pattern would suggest that size at metamorphosis 
has a large effect on fitness and that growth rates are mod- 
ulated as a correlated response. Third, developmental rates 
may increase at the expense of growth, and in this scenario 
tadpoles may speed up metamorphosis at the cost of being 
smaller at this stage. This would support a trade-off between 
growth and development, which possibly stems from a lim- 
ited resource budget that must eventually sustain either cell 
growth/division or cell differentiation (see Smith-Gill and 
Berven 1979). 

Environmental variability adds another level of complexity 
because different strategies may be favored depending on the 
conditions encountered during development. Smaller tem- 
poral ponds are expected to appear and disappear within 
a relatively short period of time, subjecting those species 
breeding in these habitats to a higher risk of desiccation. 
Evolutionary responses to increased desiccation risk may in- 
volve an overall shift toward faster development, the evolution 
of plasticity triggered only when ponds are drying (Denver 
1997), or both. Intuitively, one would expect that species from 
more temporary ponds should show an increased plasticity, 
everything else being equal, because these ponds are intrin- 
sically more variable (Newman 1992; Doughty and Reznick 
2004; Wells 2007). Conversely, the opposite pattern should 
be expected if these species are developing at rates near their 
maximum physiological capacities (Nunney 1996; Nylin and 
Gotthard 1998). It is therefore not surprising that empiri- 
cal studies testing for differences in developmental rates and 
plasticity between and within species facing contrasting des- 
iccation risks are not conclusive, and have resulted in dis- 
cordant patterns both at the interspecific and intraspecific 
level (Brady and Griffiths 2000; Leips et al. 2000; Merila et al. 
2004; Richter-Boix et al. 2006; Lind and Johansson 2007). 
Disentangling between these two alternatives involve assess- 
ing concomitantly (1) whether average developmental rates 
differ across species facing contrasting risks of desiccation 
and (2) if these species also exhibit contrasting levels of plas- 
ticity in developmental rates around these mean values. 

In the present study, we determined whether amphib- 
ian anuran species facing higher risks of desiccation have 
evolved faster developmental rates and increased plasticity 
than species living in more stable pond environments. We 
tested these hypotheses with a comparative approach and a 
large dataset that includes several anuran species belonging 
to distinct taxonomic groups. Importantly, several compar- 
ative studies have shown that evolutionary inferences may 
change when phylogenetic relations among species are taken 
into account (e.g., see Gomes et al. 2009 for a recent case 
study in anurans). We specifically examine the following 
questions: (1) Do developmental rates, mean growth rates, 
and mass at metamorphosis differ between species breeding 
in ponds with different risk of desiccation? Subsequently, we 
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determined if plastic responses showed significant interspe- 
cific differences, asking: (2) Is plasticity in developmental 
rates, mean growth rates, and mass at metamorphosis greater 
in species with increased risk of desiccation? (3) Are plas- 
tic responses in developmental rates and mean growth rates 
correlated, and do these correlations differ as a function of 
desiccation risk? Finally, we assess if evolutionary shifts in 
developmental rates have altered the magnitude and direc- 
tion of plastic responses across species: (4) Is the variation 
in mean developmental rates, mean growth rates, and body 
mass at metamorphosis observed under constant conditions 
correlated with plasticity in these variables? 

Material and Methods 
Data collection 

The literature search was performed in several electronic 
databases (ISI Web of Science, BasicBIOSIS, BioOne, Bi- 
ological Abstract, ScienceDirect, Scopus and Scirus), em- 
ploying the key words "pond desiccation," "hydroperiod," 
"metamorphosis," "life history," and "larval period," and 
the references of all papers obtained with this search were 
subsequently reviewed. Data from individual studies were 
included only if they met the following criteria. First, the 
study experimentally manipulated water level or water per- 
manency directly by comparing developmental time and size 
at metamorphosis in constant water-level treatment versus 
desiccation treatment in the laboratory or mesocosm, or in- 
directly by comparing permanent versus temporary ponds 
in the field. When more than one level was available, the 
extremes were selected to maximize the difference in desic- 
cation risk between treatments. Second, we excluded studies 
that combined pond drying effects with another variable to 
avoid confounding effects (in factorial experiments, we com- 
pared the control against the treatment where only the drying 
regime was manipulated). Third, studies that utilized snout 
vent length as an estimate of metamorphic size were not in- 
cluded, but we did include data from one study that reported 
body volume (Adams 2000) because this measurement is 
highly correlated with body mass and can be readily trans- 
formed to mass units assuming that density ~ 1 g/mL. Results 
from single studies including multiple species or populations 
were included as independent samples. 

We assembled a database comprising 62 independent cases, 
encompassing a total of 30 species, including the hybrids be- 
tween Rana sphenocephala and Rana blairi and the hybrid 
of Pelophylax lessonae and Pelophylax ridibunda (Pelophy- 
lax esculenta), that were considered as 'ecological' species 
(see Data analysis below), from 25 studies in the litera- 
ture, published between 1989 and 2006, and three unpub- 
lished datasets (Table SI). Developmental rates were calcu- 
lated as 1 /developmental period and mean growth rates as 
body mass/ developmental period. This later estimate of mean 
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growth rate was employed for simplicity as a general descrip- 
tor of growth rates (see Abrams et al. 1996, pp. 383), because 
growth trajectories in anurans can be considerably complex 
(Harris 1999) (see Fig. SI). It must be emphasized that, even 
though this estimate may be adequate for comparative pur- 
poses because it takes into consideration the final outcome of 
the tadpole's ontogeny (i.e., mass at metamorphosis), it has 
very little to say about the real shape of the developmental 
curve during this period. All subsequent statistical analyses 
have been performed employing developmental rates, though 
we often discuss results in the context of developmental peri- 
ods because it is more intuitive (e.g., Pelodytes ibericus takes 
106 days to reach metamorphic climax which corresponds to 
a developmental rate of 0.0094 days -1 , see Table SI). 

We operationally assigned species into two categories, ac- 
cording to the typical habitat they employed for breeding, 
to analyze the evolutionary consequences of developing un- 
der contrasting variation of desiccation risks: (1) High risk. 
Species exposed to a high risk of larval mortality by pond 
desiccation, including ephemeral and temporary ponds that 
hold water for only a few weeks or months and dry each year, 
occasionally with several dryings and refills per season. (2) 
Low risk. Species exposed to a lower risk of pond desiccation, 
basically permanent ponds holding water year-round in most 
years with rare events of desiccation. Desiccation risk varies 
more within and between years in ephemeral and temporary 
ponds than in permanent ponds (Richter-Boix et al. 2006). 
In most cases, we employed the habitat description provided 
in the papers to assign species/populations into these cat- 
egories. When this information was not available, we ob- 
tained species breeding habitat from AmphibiaWeb database 
(http://amphibiaweb.org/) or Global Amphibians Assess- 
ments Project database (http://www.globalamphibians.org/) 
and Lannoo (2005). We also collected information on the 
experimental venue employed in each study (i.e., laboratory, 
mesocosm or field conditions). 

Data analysis 

Analyses were performed employing phylogentic generalized 
linear models (see Garland et al. 2005 for a review), with 
the ape package available in R (http://cran.r-project.org/). 
Because we did not detect any significant effects of dif- 
ferent experimental venues in preliminary analyses, values 
obtained from laboratory, mesocosm, and field experiment 
were combined in subsequent analyses. Comparisons of de- 
velopmental strategies between High and Low risk species 
were performed employing values measured under constant 
water-level conditions, because this minimizes the poten- 
tially confounding effects of plasticity. We compared devel- 
opmental rates, mean growth, and body mass with a lin- 
ear model including desiccation risk as a categorical factor 
{High risk vs. Low risk). In these analyses, all variables were 
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logio-transformed to meet the assumption of normality. A 
similar model including developmental time as a covariable 
was subsequently employed to compare mean growth rates 
and body mass at metamorphosis controlling for differences 
in developmental time. 

To estimate and compare plastic responses across species, 
we first accounted for interspecific differences in develop- 
mental rates, mean growth rates, and body mass by setting 
these variables measured under constant conditions to 100%. 
Subsequently, we expressed estimates obtained in drying con- 
ditions as a fraction of this total, and plasticity was estimated 
as the proportional difference between mean traits under con- 
stant and drying conditions (e.g., plasticity in body mass was 
calculated as 100 x (massj — mass c ) / mass c where c and d 
correspond to constant and drying conditions, respectively). 
To ensure that none of the analyses were affected by variabil- 
ity within treatments, we also estimated plasticity as stan- 
dardized effect sizes (Hedges' d provided by meta-analysis 
performed with Metawin 2.1, Rosenberg et al. 2000), which 
was highly correlated with plasticity expressed as percentages 
(Fig. S2). Results remained qualitatively similar (analyses not 
shown), therefore we only report results with plasticity esti- 
mated as percent changes because their interpretation is more 
intuitive. Comparisons of developmental plasticity were per- 
formed with generalized linear models, as explained above, 
including desiccation risk as a categorical factor. In addi- 
tion, we ran different linear models with log 10 -transformed 
mean traits obtained in constant conditions as the indepen- 
dent variable and plasticity estimates as the dependent vari- 
able, to analyze how plasticity varied as a function of mean 
developmental and mean growth rates and body mass at 
metamorphosis. 

To conduct the phylogenetic analyses, a phylogeny at the 
family level was constructed following Frost et al. (2006). 
Importantly, the topology of this backbone phylogeny was 
congruent with more recent studies based on DNA se- 
quences (Roelants et al. 2007; Wiens 2007). Subsequently, we 
combined this information with additional detailed within- 
family phylogenetic assessments for the following families: 
Scaphiopodidae, Pelobatidae, and Pelodytidae (Garcia-Paris 
et al. 2003), Myobatrachidae (Schauble et al. 2000; Read et al. 
2001), Hylidae (Faivovich et al. 2005; Wiens et al. 2006), 
Bufonidae (Pauly et al. 2004), and Ranidae (Veith et al. 2003; 
Hillis and Wilcox 2005; Scott 2005). Multiple measurements 
per species were included as soft polytomies at the tips of the 
phylogeny, and hybrids were included in a soft polytomy with 
their respective parental species (Fig. S3). To avoid inflated 
type I error due to these polytomies (25 in total), we opted 
for a conservative approach and subtracted one degree of 
freedom for each unresolved node during hypothesis testing 
(Purvis and Garland 1993; Garland and Diaz-Uriarte 1999). 

We tested if phylogenetic signal-that is, the tendency for 
related species to resemble each other (Blomberg et al. 2003)- 
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was present for different traits as follows. We ran general- 
ized linear models employing both a star phylogeny (which 
corresponds to conventional statistical analyses, see Garland 
et al. 2005) and arbitrary branch lengths according to Pagel 
(1992), and subsequently determined how well the conven- 
tional and the phylogenetic model fitted the phenotypic data 
employing the Akaike information coefficients (AIC) (Burn- 
ham and Anderson 2002). The AIC criterion is currently a 
standard tool in model selection, which allows for comparing 
the goodness of fit of different models while penalizing for 
increasing the number of estimated parameters (the model 
with the lowest AIC value is considered the best model; Burn- 
ham and Anderson 2002). We also computed Akaike weights 
(AIC W ) to estimate the relative weight of the evidence in fa- 
vor of each model (Turkheimer et al. 2003), which can be 
loosely interpreted as the probability of each model being 
correct given all the models that were tested. These analy- 
ses indicate if signal is present in the analyzed dataset and, 
consequently, which model is more reliable for evolutionary 
inferences. 

Results 

Evolutionary differences in development 

Do developmental rates, mean growth rates, and mass at 
metamorphosis differ between species breeding in ponds 
with different risk of desiccation? Comparisons between ATC 
and AIC W show that phylogenetic models comparing devel- 
opmental rates, mean growth rates, and body mass at meta- 
morphosis had a substantially better fit than conventional 
analyses (Table 1). This suggests that these traits exhibit high 
phylogenetic signal (see Fig. 2) and those results from phylo- 
genetic analyses are more reliable. Comparisons controlling 
for phylogeny indicate that species breeding in ponds with 
increased desiccation risk have evolved significantly higher 
developmental rates than their counterparts from more per- 
manent ponds (Table 1 and Fig. 3). This difference seems to 
be primarily associated with a significantly lower body mass 
at metamorphosis in species with high risk of desiccation, 
while no statistical differences in mean growth rates were 
detected between groups (Table 1). Accordingly, when we 
control for differences in developmental rates, species from 
the High risk group exhibit significantly lower body mass at 
metamorphosis and mean growth rates (P < 0.01 in both 
cases), supporting the prediction that these species accelerate 
their development at the expense of growth (Fig. lc). 

Phenotypic plasticity 

Here we test if plasticity in developmental rates, mean growth 
rates, and mass at metamorphosis are greater in species with 
increased risk of desiccation. Comparisons between AIC val- 
ues suggest that phylogenetic models provide a better fit to 
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Table 1 . Evolutionary differences in developmental rates, mean growth rates, and body mass at metamorphosis between species breeding in ponds 
with High risk and Low risk of desiccation, measured under constant water-level conditions. Note that 25 degrees of freedom have been subtracted 
in phylogenetic analyses to account for the soft polytomies (see Methods). Models with the best fit are highlighted in bold. 



Dependent 1 


Model 


Desiccation risk 


AIC 


AIC W 2 


Developmental rates 


Conventional 


t 60 = -1.65, P = 0.104 


10.37 


0.001 




Phylogenetic 


t 35 = -2.91, P = 0.006 


-3.25 


0.999 


Mean growth rates 


Conventional 


t 60 = 0.047, P = 0.96 


114.2 


0.00 




Phylogenetic 


t 35 = 1.13, P = 0.266 


54.7 


1.00 


Mass at metamorphosis 


Conventional 


t 60 = 0.72, P = 0.474 


117.7 


0.00 




Phylogenetic 


t 35 = 3.39, P = 0.002 


37.0 


1.00 



1 Results from general linear models employing logio-transformed values (see boxplots Fig. 3). 

2 One conventional model and one phylogenetic model were analyzed per dependent variable, hence pairwise/WC„ should add up to one. 




Figure 2. Phylogenetic hypothesis employed in this study, including species/populations developmental times, body mass at metamorphosis, and 
mean growth rates under constant water level conditions. Species and population identification (see Table S1 and Fig. S3). 



plastic responses for developmental rates and body mass at 
metamorphosis, whereas conventional analyses are more re- 
liable for mean growth rates (Table S2). Results of conven- 
tional and phylogenetic models are qualitatively similar, how- 
ever, and suggest that plasticity in developmental rates, mean 
growth rates, and body mass at metamorphosis are similar 
between High and Low risk breeders (P > 0. 170 in all analyses, 
see Table S2). 



However, species breeding in ephemeral and permanent 
ponds apparently respond differently to drying conditions: 
the reduction in developmental time in species from High 
risk was accompanied by a decrease in body mass at meta- 
morphosis (Fig. 4b), whereas the largest reductions in de- 
velopment time in Low risk species were observed in species 
that increased mean growth rates (Fig. 4c). This suggests that 
a trade-off between developmental rates and mean growth 
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Figure 3. Developmental differences between species breeding in ephemeral {High risk in white) and permanent ponds (Low risk in gray), where lines 
schematically show a linear developmental trajectory with a slope corresponding to the average growth rate (see model outlined in Fig. 1 and Fig. S1). 
Boxplots depict the median, the 25% and 75% CI (box), 5% and 95% CI (error bars), and outliers for developmental period (— 1/ developmental rate), 
body mass at metamorphosis and mean growth rates obtained under constant water level conditions. Phylogenetic models always resulted in the 
best fit according to AIC values, and show that developmental rates and body mass at metamorphosis, but not mean growth rates, differ significantly 
between High risk and Low risk (Table 1 ). 



rates exists in species breeding in ephemeral ponds, while 
several species breeding in permanent ponds seem able to 
compensate accelerated developmental rates and reach meta- 
morphosis at a relatively constant body mass. These contrast- 
ing responses of Low and High risk species fit well with the 
expected scenarios proposed in Figure lb and c, respectively. 

Evolutionary shifts and their association 
with plasticity 

Does the variation in mean estimates measured under con- 
stant conditions correlate with plasticity estimates? Phyloge- 
netic generalized linear models suggest that plastic responses 
in developmental times are relatively constant across species 
and independent of the absolute length of the development 
period (f 33 = —0.62, P = 0.54). In other words, the reduction 
in developmental time due to plasticity is a relative constant 
fraction of developmental times under constant "optimal" 
conditions, regardless of whether a species takes 20dorl70d 



to develop. Conversely, the relative reduction in body mass 
at metamorphosis was significantly higher in larger species 
(f 33 = —2.86, P = 0.007, respectively). Thus, larger tadpoles 
were able to metamorphose at a substantially lower fraction 
of their mass under constant conditions. 

A similar pattern was observed for plasticity in mean 
growth rates: species growing on average faster showed larger 
reductions in mean growth rates when exposed to decreasing 
water levels. This trend was significant according to a conven- 
tional linear model (f 58 = -3.84, P = 0.0003), which fits the 
data substantially better (AIC = 537.17 and AIC W = 0.94) 
than the phylogenetic model (AIC = 543.58 and A/C w = 
0.06). Interestingly, the magnitude of the plastic response 
in developmental times was not significantly related to varia- 
tion in either body mass at metamorphosis ( f 33 = — 0.26, P = 
0.796) and mean growth rates (f 33 = 0.17, P = 0.866), hence 
the reduction in developmental times seems to be relatively 
constant across species regardless of their size at metamor- 
phosis and mean growth rates. 



20 



© 201 1 The Authors. Published by Blackwell Publishing Ltd. 



A. Richter-Boix etal. 




0 20 40 60 80 100 120 140 



Developmental time (%) 





o - 


b 








o 
irj 








!y in mas: 


o 


o 


o 


~~ o 


Plasticii 


ICC -50 
1 1 


'6' 


o ° 





-40 -20 0 

Plasticity in developmental time (%) 



c 

o 






"6" 


o °o O^^g^cgE 

o ° 


o 


i 

-40 


1 1 

-20 0 


1 

20 



Plasticity in developmental time (%) 

Figure 4. Plastic responses to pond desiccation, (a) Development dur- 
ing drying conditions compared to results observed under constant water 
level conditions that were set to 100% for comparative purposes (see 
Methods), (b) The association between plasticity in developmental time 
and plasticity in body mass at metamorphosis and (c) between plastic- 
ity in developmental time and plasticity in mean growth rate (plasticity 
estimates were calculated as 1 00 x (trait d - trait c )/ 'trait c where d and c 
correspond to drying and constant conditions, respectively). Results of 
phylogenetic regressions for species breeding in ephemeral ponds (High 
risk) and permanent ponds (Low risk) are shown with dashed and solid 
lines, respectively. The interaction between developmental time plastic- 
ity x desiccation risk was significant in analyses of mass and growth 
plasticity (P < 0.0027 in both models), indicating that species from 
High risk and Low risk employ different strategies to accelerate their 
development. 
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Discussion 

To answer the specific issues raised in the introductory sec- 
tion, our results can be succinctly summarized as follows. (1) 
Species that typically breed in temporary ponds have faster 
development rates than do species that typically breed in 
permanent ponds. This evolutionary response is associated 
with a reduced body size at metamorphosis and lower mean 
growth rates (Fig. 3). (2) There were no significant differ- 
ences in plasticity between species breeding in temporary and 
permanent ponds. Although several species were plastic and 
capable of decreasing developmental time when exposed to 
desiccation, this response was not general across taxa (Fig. 4a 
and Fig. S3). (3) Correlations between developmental rates, 
mean growth rates, and body mass at metamorphosis differed 
dramatically between species breeding in temporary and per- 
manent ponds, suggesting that evolution in developmental 
rates have resulted in changes at the level of plasticity (Fig. 4b 
and c). (4) Species with either longer larval period or larger 
body sizes at metamorphosis are not able to decrease their 
developmental period to a larger extent than other species. 

Evolutionary shifts and plastic responses 

Previous studies have suggested that anuran larvae that ex- 
hibit shorter mean developmental times in ephemeral ponds 
with increased risk of desiccation also show reduced body 
mass at metamorphosis (Denver 1997; Wells 2007). Our 
analyses show that interspecific variation in developmen- 
tal rates and body mass at metamorphosis are partly ex- 
plained by differences in desiccation risk during larval de- 
velopment (Table 1, Fig. 3). Species breeding in ephemeral 
ponds enter metamorphosis at a comparatively smaller body 
size than their counterparts breeding in permanent ponds, 
which seems to explain the significant difference in devel- 
opmental period observed between groups. If a critical size 
is necessary to trigger metamorphosis, as proposed by some 
theoretical models (Wilbur and Collins 1973; Day and Rowe 
2002), this would suggest that High risk species have evolved 
lower threshold sizes to reduce developmental time. Morey 
and Reznick (2004) reported that spadefoot toad species 
inhabiting ephemeral ponds had the shortest development 
times and threshold size to enter metamorphosis, and our 
results suggest that this is a general pattern. 

In addition, analyses controlling for differences in develop- 
mental rates suggest that species breeding in ephemeral ponds 
actually grow slower than their counterparts breeding in per- 
manent ponds. This is counterintuitive if one assumes that 
a critical size is necessary to trigger metamorphosis into the 
adult form, given that lower overall growth rates would delay 
metamorphosis and may potentially have a negative impact in 
fitness (Day and Rowe 2002). Taken together, developmental 
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differences observed between groups indicate that increased 
pressures on larval survival in ephemeral ponds may have 
resulted in accelerated developmental rates at the expense of 
growth, which is adaptive if the benefits of developing faster 
and evolving lower threshold sizes are higher than the costs 
they impinge on growth rates. Importantly, these differences 
were detected under constant water-level conditions, suggest- 
ing that evolutionary responses to increased desiccation risks 
partly involve an evolutionary shift in mean values toward 
more accelerated development. 

At the level of developmental plasticity, comparisons be- 
tween species breeding in ephemeral and permanent ponds 
showed that both groups showed a similar reduction in de- 
velopmental period in response to drying conditions. How- 
ever, the nature of plastic responses actually depended on 
the breeding habitat: whereas species breeding in ephemeral 
ponds showed pronounced reductions in size at metamor- 
phosis in response to drying conditions, this was generally not 
observed in species inhabiting permanent ponds (Fig. 4b). 
Even though a reduction in mean growth rate may be partly 
explained by increased stressful conditions when ponds are 
drying, this should not result in significant differences be- 
tween species breeding in ephemeral and permanent ponds. 
Instead, contrasting patterns between groups suggest that 
species breeding in ephemeral ponds allocate resources pref- 
erentially to development rather than growth (Figs. 1 c and 4c; 
Wilbur and Collins 1973; De Witt et al. 1998; Harris 1999). 



Evolutionary trade-offs and constraints 

Development and growth involve two very distinct physio- 
logical processes, where the former is essentially associated 
with cell differentiation while the second is primarily de- 
termined by cell proliferation and growth (Smith-Gill and 
Berven 1979). Although these processes are intrinsically con- 
nected, some degree of independence between them must 
ultimately account for the diversity of developmental strate- 
gies observed across anuran species (Fig. 3), as predicted by 
Wilbur and Collin's (1973) model. Our analyses are partic- 
ularly relevant in this context because species under strong 
selection for faster development during larval stages appar- 
ently maximize their developmental rates at the expense of 
growth, which suggests a trade-off between development and 
growth (Fig. lc). Several lines of evidence support this con- 
clusion. For instance, these species exhibit significantly slower 
mean growth rates in analyses controlling for differences in 
developmental rates under constant water-level conditions. 
In addition, the association between developmental rates and 
mean growth rates is nearly flat in this group, contrasting with 
the positive association between these variables observed in 
species breeding in permanent ponds, as indicated by the 
nearly significant interaction in the general linear model. 



Furthermore, contrasting plastic responses to drying con- 
ditions indicate that species breeding in ephemeral ponds 
cannot increase mean growth rates under stressful condi- 
tions, which is expected if mean growth rates under constant 
water-level conditions are close to a physiological limit. The 
maximum observed increase in mean growth rates due to 
plasticity was 6.9% in species from ephemeral ponds and 
74.5% in species breeding in permanent ponds. The mean 
response within the sample that showed an increase in mean 
growth rates ( n = 16 cases) provided a similar picture, an av- 
erage (±SE) increase of 3.6 ± 0.6% in species from ephemeral 
ponds versus a 29.1 + 8.3% increase in their counterparts 
from permanent ponds (/- 14 = 3.5, P = 0.0036). This pro- 
vides compelling evidence that anuran species breeding in 
ephemeral ponds are maximizing mean growth rates and con- 
sequently show little plasticity in this trait (Newman 1 988; Re- 
ques and Tejedo 1997). Conversely, long-lasting pond breed- 
ers can increase mean growth rates under desiccation risks, 
hence mean growth rates under constant water-level condi- 
tions are apparently not being maximized (see also Wilbur 
1987; Semlitsch and Wilbur 1988; Boone et al. 2004). Inter- 
estingly, some clades with facultative or obligate carnivory 
within the families Scaphiopodidae (Pfenning 1992) and 
Ceratophryidae (Cei 1980) breed in ephemeral ponds and 
have very short larval periods (Buchholz and Hayes 2002; 
Marangoni et al. 2009), suggesting that these diet shifts might 
have resulted from selection on increased mean growth rates. 

While most species breeding in ephemeral ponds showed 
similar responses to drying conditions (a reduction in both 
mass at metamorphosis and mean growth rate), considerable 
variation in the nature and the magnitude of plastic responses 
was observed across species breeding in permanent ponds. 
Many species did not increase their developmental rates in 
response to drying conditions, even though some of these 
species exhibited a reduction in body size at metamorphosis 
of almost 40% (Fig. 4b) . Conversely, other species showed im- 
portant reductions in developmental time and were also able 
to increase mean growth rates as a compensatory response 
(Fig. 4c), suggesting that the primary target of selection in 
these species is body mass at metamorphosis. Importantly, 
this response was observed in distantly related taxa (Fig. S3), 
which explains why conventional statistics assuming a star 
phylogeny, provided a better fit in analyses of mean growth 
rates. Patterns observed in this subset agree with observa- 
tions carried out in butterflies and damselfiies, which exhibit 
a shortening in developmental period while simultaneously 
reducing size and accelerating mean growth rates (Nylin et al. 
1996; Strobbe and Stoks 2004; De Block et al. 2008). In 
summary, our results show that plastic responses are con- 
strained in species with accelerated developmental rates and 
suggest that they are indeed developing near their maximum 
physiological capacities. A similar comparative approach may 
shed light on the generality of this observation across other 
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biological systems in which developmental rates are expected 
to be under selection. 

Phylogenetic effects 

Absolute limits to performance certainly exist, but this does 
not account for the interspecific differences observed in de- 
velopmental strategies. Although differences may be partly 
due to the stressful physical and biotic conditions of desiccat- 
ing ponds, our results suggest that an important fraction of 
the variation can be attributed to phylogenetic history. For ex- 
ample, bufonids in general metamorphose at very small body 
sizes, whereas the opposite pattern is true for scaphiopodids 
(Fig. 2, see also Werner 1986). It is possible that contrast- 
ing larval sizes reflect differences in fecundity across species, 
hence more detailed information on maternal investment, 
egg, and clutch size may shed light on the relevant mecha- 
nisms underlying developmental differences between species 
(Wells 2007, pp. 494-515). In addition, physiological con- 
straints may differ across taxa (e.g., rudimentary nonfunc- 
tional lungs are characteristic of bufonid larvae; Ultsch et al. 
1999 and references therein). 

Understanding which factors ultimately explain these dif- 
ferences is crucial to determine if and how anuran species 
may respond to increasing risks of desiccation. Even though 
comparative studies are strictly correlational and may provide 
limited information on the mechanisms underlying species 
developmental differences, phylogenetic information can be 
valuable for predictive purposes, because closely related taxa 
seem to employ similar developmental strategies and may 
also potentially share the same physiological limitations. 

Acknowledgments 

We are very grateful to M. D. Boone, J. Loman for provid- 
ing important unpublished data, F. Johansson, A. Laurila, C. 
Navas, R. Stoks, J. Wiens, and two anonymous reviewers for 
their comments on earlier versions of the manuscript. This re- 
search was partially funded by project (CGL2004-0 1872/BOS 
and CGL2009-12767-C02-02) from Direction General de In- 
vestigation Cientifica y Tecnica (MICINN), conceded to M. 
Tejedo and by the grant BFU2009-07564 from the Ministe- 
rio de Ciencia e Innovation (MICINN, Spain) awarded to 
ELR. AR-B was supported by a Spanish Ministry of Educa- 
tion and Culture postdoctoral grant (MEC2007-0944). ELR 
is currently a Ramon y Cajal Fellow funded by the Ministerio 
de Ciencia e Innovation (Spain). 

References 

Abrams, P. A., O. Leimar, S. Nylin and C. Wiklund. 1996. The 
effect of flexible growth rates on optimal sizes and 
development times in a seasonal environment. Am. Nat. 
147:381-395. 



Evolution and Plasticity of Anuran Development 

Adams, M. J. 2000. Pond permanence and the effects 

of exotic vertebrates on anurans. Ecological Applications, 
10:559-568. 

Arendt, J. D. 1997. Adaptive intrinsic growth rates: an integration 

across taxa. Quat. Rev. Biol. 72:149-177. 
Babbitt, K. J., M. J. Baber and T. L. Tarr. 2003. Patterns of larval 

amphibian distribution along a wetland hydroperiod gradient. 

Can. J. Zool. 81:1539-1552. 
Blomberg, S. P., T. Garland, Jr. and A. R. Ives. 2003. Testing for 

phylogenetic signal in comparative data: behavioral traits are 

more labile. Evolution 57:717-745. 
Boone, M. D., E. E. Little and R. D. Semlitsch. 2004. Overwintered 

bullfrog tadpoles negatively affect salamanders and anurans in 

native amphibian communities. Copeia 2004:683-690. 
Brady, L. D. and R. A. Griffiths. 2000. Developmental responses 

to pond desiccation in tadpoles of the British anuran 

amphibians {Bufo bufo, B. calamita and Rana temporaria). 

J. Zool. 252:61-69. 
Buchholz, D. R. and T. B. Hayes. 2002. Evolutionary patterns of 

diversity in spadefoot toad metamorphosis (Anura: 

Pelobatidae). Copeia 2002:180-189. 
Burnham, K. P. and D. R. Anderson. 20 02. Model selection and 

multi-model inference: a practical information-theoretic 

approach. 2nd edition. Springer, New York. 
Cei, J. M. 1980. Amphibians of Argentina. Monitore Zoologico 

Italiana (N.S.). 243 Monografia 2:1-609. 
Day, T. and L. Rowe. 2002. Developmental thresholds and the 

evolution of reaction norms for age and size at life-history 

transitions. Am. Nat. 159:338-350. 
De Block, M., M. A. McPeek and R. Stoks. 2008 . Life-history 

evolution when Lestes damselflies invaded vernal ponds. 

Evolution 62:485-493. 
Denver, R. J. 1997. Proximate mechanisms of phenotypic 

plasticity in amphibian metamorphosis. Am. Zool. 

37:172-184. 

De Witt, T. J., A. Sih and D. S. Wilson. 1998. Costs and limits of 
phenotypic plasticity. Trends Ecol. Evol. 13:77-81. 

Doughty, P. and D. N. Reznick. 2004. Patterns and analysis of 
adaptive phenotypic plasticity in animals. In: T. DeWitt and S. 
M. Scheiner, editors. Phenotypic plasticity: functional and 
conceptual approaches. Oxford University Press, Oxford, 
pp. 126-150. 

Faivovich, J., C. F. B. Haddad, P. C. A. Garcia, D. R. Frost, J. A. 

Campbell and W. C. Wheeler. 2005. Systematic review of the 

frog family Hylidae, with special reference to Hylinae: 

phylogenetic analysis and taxonomic revision. Bull. Am. Mus. 

Nat. Hist. 294:1-240. 
Frost, D. R., T. Grant, J. Faivovich, R. H. Bain, A. Haas, C. F. B. 

Haddad, R. O. de Si, A. Channing, M. Wilkinson, S. C. 

Donnellan. 2006. The amphibian tree of life. Bull. Am. Mus. 

Nat. Hist. 297:1-371. 
Garcia-Paris, M., D. R. Buchholz and G. Parra-Olea. 2003. 

Phylogenetic relationships of Pelobatoidea re-examined using 

mtDNA. Mol. Phylogenet. Evol. 28:12-23. 



© 201 1 The Authors. Published by Blackwell Publishing Ltd. 



23 



Evolution and Plasticity of Anuran Development 

Garland, T., Jr., A. F. Bennett and E. L. Rezende. 2005. 

Phylogenetic approaches in comparative physiology. J. Exp. 

Biol. 208:3015-3035. 
Garland, T. Jr., and R. Diaz-Uriarte. (1999). Polytomies and 

independent contrasts: an examination of the bounded 

degrees of freedom approach. Systematic Biology, 48:547-558. 
Gomes, F. R., E. L. Rezende, M. B. Grizante and C. A. Navas. 

2009. The evolution of jumping performance in anurans: 

morphological correlates and ecological implications. J. Evol. 

Biol. 22:1088-1097. 
Gotthard, K. and S. Nylin. 1995. Adaptive plasticity and plasticity 

as an adaptation: a selective review of plasticity in animal 

morphology and life history. Oikos 75:3-17. 
Harris, R. N. 1999. The anuran tadpole: evolution and 

maintenance. In: R. W. McDiarmid and R. Altig, editors. 

Tadpoles: the biology of anuran larvae. University of Chicago 

Press, Chicago, pp. 279-294. 
Hillis, D. M. and T. P. Wilcox. 2005. Phylogeny of the New World 

true frogs {Rand). Mol. Phylogenet. Evol. 34:299-314. 
Johansson, E, R. Stoks, L. Rowe and M. De Block. 2001. Life 

history plasticity in a damselfly: effects of combined time and 

biotic constraints. Ecology 82:1857-1869. 
Lannoo, M. 2005. Amphibian declines: the conservation status of 

United States species. University of California Press, Berkeley, 

California, USA. 
Leips, J., M. G. McManus and J. Travis. 2000. Response of 

treefrog larvae to drying ponds: comparing temporary and 

permanent pond breeders. Ecology 81:2997-3008. 
Lind, M. I. and F. Johansson. 2007. The degree of adaptive 

phenotypic plasticity is correlated with the spatial 

environmental heterogeneity experienced by island 

populations of Rana temporaria. J. Evol. Biol. 20:1288-1297. 
Marangoni, E, E. Schaefer, R. Cajade and M. Tejedo. 2009. 

Growth marks formation and chronology of two neotropical 

anuran species. J. Herpetol. 43:546-550. 
Merila, J., A. Laurila and B. Lindgren. 2004. Variation in the 

degree and costs of adaptive phenotypic plasticity among Rana 

temporaria populations. J. Evol. Biol. 17:1132-1140. 
Morey, S. and D. N. Reznick. 2004. The relationship between 

habitat permanence and larval development in California 

spadefoot toads: field and laboratory comparisons of 

developmental plasticity. Oikos 104:172-190. 
Newman, R. A. 1988. Adaptive plasticity in development of 

Scaphiopus couchii tadpoles in desert ponds. Evolution 

42:774-783. 

Newman, R. A. 1992. Adaptive plasticity in amphibian 

metamorphosis. Bioscience 42:671-678. 
Nunney, L. 1996. The response to selection for fast larval 

development in Drosophila melanogaster and its effect on adult 

weight: an example of a fitness trade-off. Evolution 

50:1193-204. 

Nylin, S., K. Gotthard and C. Wiklund. 1996. Reaction norms 
and size at maturity in Lasiommata butterflies: predictions and 
tests. Evolution 50:1351-1358. 



A. Richter-Boix etal. 

Nylin, S. and K. Gotthard. 1998. Plasticity in life-history traits. 

Annu. Rev. Entomol. 43:63-83. 
Pagel, M. 1992. A method for the analysis of comparative data. 

J. Theor. Biol. 156:431-442. 
Pauly, G. B., D. M. Hillis and D. C. Cannatella. 2004. The history 

of a nearctic colonization: molecular phylogenetics and 

biogeography of the nearctic toads (Bufo). Evolution 

58:2517-2535. 

Pfenning, D. W. 1992. Proximate and functional causes of 
polyphenism in an anuran tadpole. Funct. Ecol. 6:167-174. 

Purvis, A., and T. Garland Jr. (1993. Polytomies in comparative 
analyses of continuous data. Systematic Biology, 
42:569-575. 

Read, K., J. Scott Keogh, I. A. W. Scott, J. Dale Roberts and P. 

Doughty. 2001. Molecular phylogeny of the Australian frog 

genera Crinia, Geocrinia, and allied taxa (Anura: 

Myobatrachidae). Mol. Phylogenet. Evol. 21:294-308. 
Reques, R. and M. Tejedo. 1997. Reaction norms for 

metamorphic traits in natterjack toads to larval density and 

pond duration. J. Evol. Biol. 10:829-851. 
Richardson, J. M. L. 2001. The relative roles of adaptation and 

phylogeny in determination of larval traits in diversifying 

anuran lineages. Am. Nat. 157:282-299. 
Richter-Boix, A., G. A. Llorente and A. Montori. 2006. A 

comparative analysis of the adaptive developmental hypothesis 

in six Mediterranean anuran species along a pond permanency 

gradient. Evol. Ecol. Res. 8:1139-1154. 
Roelants, K., D. J. Gower, M. Wilkinson, S. P. Loader, S. D. Biju, 

K. Guillaume, M. Linde and F. Bossuyt. 2007. Global patterns 

of diversification in the history of modern amphibians. Proc. 

Natl. Acad. Sci. U. S. A. 104:887-892. 
Rosenberg, M. S., D. C. Adams and J. Gurevitch. 2000. MetaWin: 

statistical software for meta-analysis. Version 2.1. Sinauer 

Associates, Inc., Sunderland, Massachussetts, USA. 
Schauble, C. S., C. Moritz and R. W. Slade. 2000. A molecular 

phylogeny for the frog genus Limnodynastes (Anura: 

Myobatrachidae). Mol. Phylogenet. Evol. 16:379-391. 
Scott, E. 2005. A phylogeny of ranid frogs (Anura: Ranoidea: 

Ranidae) based on a simultaneous analysis of morphological 

and molecular data. Cladistics 21:507-574. 
Semlitsch, R. D. and H. M. Wilbur. 1988. Effects of pond drying 

time on metamorphosis and survival in the salamander 

Ambystoma talpoideum. Copeia 1988:978-983. 
Smith-Gill, S. J. and K. A. Berven. 1979. Predicting amphibian 

metamorphosis. Am. Nat. 113:563-585. 
Stearns, S. C. and J. C. Koella. 1986. The evolution of phenotypic 

plasticity in life-history traits: predictions of reaction norms 

for age and size at maturity. Evolution 40:893-913. 
Stoks, R., M. De Block and M. A. McPeek. 2006. Physiological 

costs of compensatory growth in a damselfly. Ecology 

87:1566-1574. 

Strobbe, F. and R. Stoks. 2004. Life history reaction norms to 
time constraints in a damselfly: differential effects on size and 
mass. Biol. J. Linn. Soc. Lond. B 83:187-196. 



24 



© 201 1 The Authors. Published by Blackwell Publishing Ltd. 



A. Richter-Boix etal. 



Evolution and Plasticity of Anuran Development 



Turkheimer, F. E., R. Hinz and V. J. Cunningham. 2003. On the 
undecidability among kinetic models: from model selection to 
model averaging. J. Cereb. Blood Flow Metab. 
23:490-498. 

Ultsch, G. R., D. F. Bradford, and J. Freda. 1999. Physiology: 

coping with the environment. In: R. W. McDiarmid and R. 

Altig, editors. Tadpoles: the biology of anuran larvae. 

University of Chicago Press.Chicago. pp 189-214. 
Veith, M., J. Kosuch, and M. Vences. 2003. Climatic oscillations 

triggered post-Messinian speciation of Western Paleartic 

brown frogs (Amphibia, Anura, Ranidae). Phylogenetics and 

Evolution, 26:310-327. 
Wells, K. D. 2007. The ecology and behavior of amphibians. 

University of Chicago Press, Chicago. 
Werner, E. E. 1986. Amphibian metamorphosis: growth rate, 

predation risk, and the optimal size at transformation. Am. 

Nat. 128:319-341. 
Wiens, J. J., C. H. Graham, D. S. Moen, S. A. Smith and T. W. 

Reeder. 2006. Evolutionary and ecological causes of the 

latitudinal diversity gradient in Hylid frogs: treefrog trees 

unearth the roots of high tropical diversity. Am. Nat. 

168:579-596. 

Wiens, J. J. 2007. Global patterns of diversification and species 
richness in amphibians. Am. Nat. 170:S86-S106. 



Wilbur, H. M. 1987. Regulation of structure in complex 

systems: experimental temporary pond communities. 

Ecology 68:1437-1452. 
Wilbur, H. M. and J. P. Collins. 1973. Ecological aspects of 

amphibian metamorphosis. Science 182:1305-1314. 
Woodward, B. D. 1983. Predator-prey interactions and 

breeding-pond use of temporary-pond species in a desert 

anuran community. Ecology 64:1549-1555. 

Supporting Information 

Additional Supporting Information may be found online on 
Wiley Online Library. 

Table SI. Dataset 

Table S2. Regression models. 

Figure SI. Growth trajectory of two hypothetical species and 

definition of mean growth rate. 

Figure S2. Developmental plasticity estimates. 

Figure S3. Phylogenetic hypothesis and phenotypic plasticity. 

Please note: Wiley- Blackwell is not responsible for the content 
or functionality of any supporting materials supplied by the 
authors. Any queries (other than missing material) should be 
directed to the corresponding author for the article. 



© 201 1 The Authors. Published by Blackwell Publishing Ltd. 



25 



